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A microscopic calculation of reaction cross sections for nucleon-nucleus scattering has been per- 
formed by explicitly coupling the elastic channel to all particle-hole excitations in the target and 
one-nucleon pickup channels. The particle-hole states may be regarded as doorway states through 
which the flux flows to more complicated configurations, and subsequently to long-lived compound 
nucleus resonances. Target excitations for '*°''**Ca, ^*Ni, ^"Zr and ^^''Sm were described in a random- 
phase framework using a Skyrme functional. Reaction cross sections obtained agree very well with 
experimental data and predictions of a state-of-the-art fitted optical potential. Couplings between 
inelastic states were found to be negligible, while the pickup channels contribute significantly. The 
effect of resonances from higher-order channels was assessed. Elastic angular distributions were also 
calculated within the same method, achieving good agreement with experimental data. For the first 
time observed absorptions are completely accounted for by explicit channel coupling, for incident 
energies between 10 and 70 MeV, with consistent angular distribution results. 

PACS numbers: 24.10.-i, 24.10.Eq, 24.50.-|-g, 25.40.-h 



I. INTRODUCTION 

Achieving a quantitative and predictive description of 
the structure of and reactions wiili nuclei across the iso- 
topic chart is an important and challenging goal of nu- 
clear physics. Accurate knowledge of reaction rates, in 
particular those related to reactions induced by a single 
nucleon, are important for nuclear energy applications 
and for understanding astrophysical phenomena pLj , such 
as the evolution of stars and the synthesis of the elements. 
Radiobiology and space science developments rely on the 
proper determination of reaction observables to provide 
accurate yields and spectra for radiation protection and 
risk estimates [5]. National security applications also 
make use of reaction and structure information to detect 
nuclear materials of interest. 

Accurate prediction of quantities related to nuclear re- 
actions is a complex problem as not only the desired out- 
come, i.e. exit channel, has to be considered, but also 
the interference and competition with all other possi- 
ble outgoing channels. A successful account of elastic 
nucleon-nucleus scattering, for example, has to include 
the effects from the excitation of non-elastic degrees of 
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freedom, such as collective and particle-hole (p-h) exci- 
tations, transfer reactions, etc. The picture that emerges 
is one in which flux is removed from the elastic channel 
by couplings to the non-elastic degrees of freedom. For- 
mally, these non-elastic effects can be accounted for by 
the projection-operator approach of Feshbach [5]. This 
approach reduces the complexity of the problem by in- 
troducing an effective optical model potential (OMP). 
This potential, often complex, can be defined [31 3] as the 
effective interaction in a single-channel calculation that 
contains the effects of all the other processes that occur 
during collisions between nuclei. OMPs play a very im- 
portant role in the description of nuclear reactions. They 
are extensively used to describe the interactions of pro- 
jectile and target in the entrance channel, and the inter- 
action of ejectile and residual nuclei after the reaction; 
they are crucial ingredients in direct-reaction analyses 
(e.g. elastic and inelastic scattering, transfer reactions, 
etc.) and provide transmission coefficients for statistical 
(Hauser-Feshbach) calculations. 

Most widely used are phenomenological OMPs fitted 
to reproduce experimental data sets. They have been 
extremely successful for many applications involving nu- 
clei in the range of the fits [5]. At the same time, such 
adjustable potentials make strong assumptions about lo- 
cality and momentum dependence that are probably not 
justified. In addition, present nuclear theory applica- 
tions require increasingly accurate predictions, specially 
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for isotopes off stability. However, for nuclei lying outside 
the range of the fits, such as unstable nuclei produced at 
rare-isotope facilities, in the r-process, and in advanced 
reactor applications, this can lead to unquantifiable un- 
certainties. The reaction mechanisms of systems involv- 
ing weakly-bound nuclei are known to be different from 
the strongly-bound ones [HI [Zj . Studies have shown that 
the behaviour of phenomenological OMPs for stable sys- 
tems usually cannot be directly extended to systems in- 
volving weakly-bound nuclei [8] . To achieve a better un- 
derstanding of nuclear reactions and structure, it is thus 
important to calculate OMPs by more fundamental, gen- 
eral, and first-principle methods. 

According to microscopic reaction theory, an OMP is 
comprised of two components. The first is a real bare 
potential, corresponding to the diagonal potential within 
the elastic channel, which is generally obtained by folding 
the nucleon distributions of both nuclei with a nucleon- 
nucleon effective interaction. The second is a complex dy- 
namic polarization potential which arises from couplings 
to inelastic states. The resulting optical potential is com- 
plex, composed of a real part (usually slightly different 
from the bare potential) and an imaginary component. 
The latter gives rise to absorption of flux from the elas- 
tic channel to the other reaction channels, and is hence 
directly connected with observed reaction cross-sections. 

Several attempts have been made to generate OMPs 
from microscopic approaches. Some have used the nu- 
clear matter approach j^, where the calculation is first 
performed in nuclear matter and a potential is then ob- 
tained for finite nuclei by an appropriate local density 
approximation. This approach provides accurate results 
at nucleon energies > 50 MeV [T^. Recently, new meth- 
ods based on self-energy theory have been implemented 
[llj . and new calculations, which combine a nuclear mat- 
ter approach and Hartree-Fock-Bogoliubov (HFB) mean 
field structure model, provide encouraging results for 
neutron scattering below 15 MeV [TH |T3]. Earlier at- 
tempts used the nuclear structure approach, which is 
more suitable at energies below 50 MeV [M], and calcu- 
lated second-order diagrams using particle-hole propaga- 
tors in the random-phase approximation (RPA) [T^17| . 
However, these were not able to fully explain observed 
absorption: e.g., in Ref. [16], the couplings could account 
only for « 44% of the nucleon-nucleus absorption and, in 
Ref. [T7], only for « 71% including charge exchange. 

The construction of OMPs from microscopic methods 
that use mapping of effective interactions to nucleon- 
nucleon g matrices, which are solutions of nuclear mat- 
ter equations, have proven to provide good agreement 
with reaction cross-section data [THlIin]. However, due to 
the increased number of discrete states of heavier nuclei 
[A > 12), this approach does not accommodate their spe- 
cific structure and processes. A method that implements 
successive spectator expansions of the optical potential, 
where the projectile is considered to interact at first or- 
der with a single nucleon of the target, has also proven 
to be a successful tool to obtain nucleon-nucleus cross 



sections [20] i although analyses were made only for scat- 
tering energies in the range ~65 to 400 MeV. In addition, 
the method described in Ref. [5D] treats the propagator 
modification through a nuclear mean field potential taken 
from structure calculations. Although this is a valid ap- 
proach, it may not be completely satisfactory when aim- 
ing to keep full consistency with the theory of multiple 
scattering. A more consistent alternative, although not 
intractable, would be more difficult to implement. 

RPA-based microscopic methods have also been ap- 
plied to heavy-ion reactions achieving good description of 
the double giant dipole resonance [2Tti23] . In such works, 
one- and two-phonon states were populated and anhar- 
monicities and nonlinear terms were treated. OMPs were 
then obtained through a semiclassical approach, by inte- 
grating the excitation probability over all impact param- 
eters. 

The latest advancements in the description of the 
structure of the nuclei from ah initio methods [24 allow 
for the development of more fundamental reaction mod- 
els based directly on structure results. Such microscopic 
approach leads to the calculation of reaction observables 
that are consistent with the structure inputs adopted. 
In addition, OMPs based on microscopic approaches are 
much more reliable than phenomenological ones when ex- 
trapolated to describe processes involving unstable nu- 
clei or previously unquantified transitions. Among the 
different microscopic structure models, methods based 
on energy-density functional (EDF) emerge as the only 
tractable theoretical tool that can be applied to all the 
nuclides with ^ > 40 [24j . 

In this paper we report on recent progress made to- 
wards achieving a complete microscopic calculation of 
the reaction cross sections for both neutron and proton 
induced reactions on a variety of medium-mass targets. 
We make extensive use of recently developed fundamen- 
tal structure models based on energy-density-functional 
theory, such as RPA and quasi-particle RPA (QRPA). 
We calculate sets of excited states and the corresponding 
transition densities and potentials from the ground state 
(g.s.) for different nuclei across the periodic table. We 
then incorporate this information about the transitions in 
coupled-reaction channels calculations of nucleon-nucleus 
reactions, coupling to all relevant channels necessary to 
consistently obtain accurate cross sections. First results 
of this approach were reported in Ref. 525] , while here we 
detail and extend that work. 

The paper is structured as follows: in Section jll] we 
present details of the structure models used and explain 
how they connect with reaction calculations; Section jllT] 
explains the procedure adopted for the coupled reaction 
channel calculations; in Section |IV| our main results are 
shown and discussed; and Section |V] presents our final 
conclusions. 
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II. STRUCTURE MODELS 

A HFB calculation gives the particle and hole levels 
of a given nucleus and fixes the p-h basis states for gen- 
erating excited states within the framework of (Q)RPA, 
thus accounting for long range correlations caused by the 
residual interactions within the target. To obtain the ini- 
tially occupied proton and neutron levels in a nucleus, we 
use the Skryme energy-density functional SLy4 [211 Table 
1], a parametrization designed to describe systems with 
arbitrary neutron excess, from stable to neutron matter, 
by improving isotopic properties, which overcomes defi- 
ciencies of other interactions away from the stability line. 
Although we used only the SLy4 parametrization in our 
work, the method is general enough to use any Skyrme 
force or any other functional. 

A. Ingredients of the Calculation 

In this subsection, we show the equation of the single- 
particle wave function and building block of the RPA 
excited state. To describe the Hartree-Fock basis formed 
by particles with spin s = |, the following state vectors 
can be defined: 

Inl^jnit) = ^' [Yi{r) x xihm xY°, (1) 

with radial wave functions (j)nijt{f) expressed in a co- 
ordinate representation, and xi and xY° are spin and 

isospin components, respectively. We associate the cre- 
ation and annihilation operators a}^igj^i and anisjmt with 
these state vectors. It is convenient to define a modified 
annihilation operator, anisjmt = (-1)-'""' anZsjmi, which 
is a spherical tensor of rank j and projection —to. 

We now define particle {p) and hole [h) states corre- 
sponding to orbitals above or below the Fermi surface, 
schematically indicated hy p > F and h < F, respec- 
tively. The symbol p represents all quantum numbers 
except the magnetic projection; i.e. p = {uplp^jptp}. 
The same definition applies to the hole states, replacing 
p by h. 

With the above definitions we define an operator that 
creates a particle-hole pair coupled to angular momen- 
tum J and projection A/, 

(P' ^ i^pmpjh ~mh\J M) al„,^ Uhm,, . (2) 

The Hermitian conjugate of h), 

Ajm{p, h) - (-1)^-*' Aj-m{p, h), (3) 

destroys a particle-hole pair, where M = — M. Both 
A^jf^jijp, h) and Ajf^jip, h) are spherical tensors of rank J 
and projection M. 

The QRPA formalism extends this HFB approach by 
including quasiparticle excitations in the following way: 



1) the single-particle wave function is extended to a two- 
component representation [27] ; 2) the particle-creation- 
hole-annhilation and particle-annihilation-hole-creation 
are extended to two-quasiparticle creation and annihi- 
lation ^28n29j. 

The HFB calculations were performed using a slightly 
modified version of the hfbrad ^ code called hfb- 
MARIO (version 6.2) |30j. 



B. RPA States 

We define a boson operator Ojyjj^^ that creates the 
RPA state \NJM), where TV is the principal quantum 
number, when applied to the correlated ground state |0), 
which is the vacuum for the RPA excitations. We also 
define the time-reversed destruction operator 

0^JM = (-1)''^*' Qnj-m. (4) 

The particle-hole operators relate to the boson operators 
through 

^NJM= E A^AliPM-Y^'',' AjM{p,h), 

p>F, h<F 

(5) 

^NJM — E '^pi^^ Ajm{p, h) — Yj^'' A^j]^,j{p, h), 

p>F. h<F 

(6) 

where X^ff and Y^'^ are the components of the linear 
combination of particle- hole excitations used to construct 
the RPA states [28l [31] . In the present treatment the X 
and Y coefRcients are real and obey the normalization 
and orthogonality relations 

^ph ^ph - ^ph J^ph —ONN' [i) 

p>F, h<F 

and 

^ph "ph ~ ^ph ^ph — \°J 

p>F, h<F 

Such X and Y amplitudes were obtained by solving the 
RPA equations [32] with the particle-hole interaction de- 
rived from the second derivative of the energy density 
functional 1331 . 



C. RPA Transition Densities 

Transitions of the many-nucleon system from an initial 
RPA state \i) = \aiIiMi), with angular momentum 
and projection Mi to a final RPA state |/) = \afIfMf), 
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where and a/ are additional quantum numbers re- 
quired to characterize the states, can be expressed with 
the help of transition densities [31] : 

Tq, flf N 

Psu y^t) 



An 



{Lf,S,^\JM) {UM,JM\IfMf)pl'i^/ (r)Yl^ (f , 

(9) 

where the sum over index n represents summation over 
all occupied orbitals (for protons or neutrons), is the 
position at which the transition density is calculated, and 
r„ is the spatial-coordinate operator of the n'th particle 
in the Hilbert space containing the states. The Sg^, and 
TrPg are the spin and isospin transition operators respec- 
tively; T = corresponds to the isoscalar part of the 
interaction and T = 1 to the isovector one. For 5 = 
or 1, 

Sqo = 1 and Sii^ = ui^, (10) 
together with similar quantities Trg for isospin, 

7^)0 = 1 and Tlq = T-yq, (11) 

where a\i, and T\q are the spherical components of the 
vector of Pauli matrices. It is convenient to introduce 
a multipole expansion for the transition densities. The 
coordinate-space radial multipole transition density is 

(12) 

Using second-quantization techniques, we can write Eq. 
as 



2L 



1 



2// + 1 E 

aijltl 



X {a2j2t2\\fL{r)[YL{r) x Ss]j\\aditi) , 

(13) 

where the isospin matrix element is 7^^*^ ~ 
{■^t2\TTq\-^ti) . In this expression the density is a sum over 
reduced matrix elements between single-particle states 
weighted by spectroscopic amplitudes Z^^^^{a2j2,Oiiji), 
the calculation of which will be described below. 

We can calculate the corresponding momentum-space 
density by the Fourier Bessel transform, as described in 
Appendix [Xj obtaining 

(14) 

X (a///ll E ^Liqrn) [^^(fn) x S^]j T^JIa./,). 



Transition densities can be classified according to the 
specific states that they connect. For a spin-0 even-even 
spherical nucleus, we can divide the problem of finding 
such densities into three parts, according to the number 
of RPA phonons in the initial and final states: 1. No 
phonons in either initial or final state (this is needed 
for elastic scattering from the ground state, described in 
Sec. |II C~ ); 2. No phonon in the initial state, one phonon 
in the final state (this is needed for inelastic scattering 
from the ground state to an excited state, described in 
Sec. II C 2 1 ; 3. One phonon in each of the initial and final 



states (this corresponds to inelastic scattering between 
two excited states, and also to elastic scattering from an 
excited state. 
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It is described in Sec. 
In each case the main work is the calculation of the 
spectroscopic amplitudes (or Z coefficients). Using these 
we construct the radial densities needed for folding-model 
calculations of transition and diagonal potentials. Of 
course, the same coefficients can be used to calculate any 
other one-body operator connecting the initial and final 
states. 

We use the symbol p to represent the set of quantum 
numbers {riplp^jptp} for particle states, along with an 
equivalent definition for hole states. When the particle 
has a subscript (e.g. pi), rather than use the cumber- 
some pi = ^jpiipil we use the simpler notation 
{nili^jiti} unless it leads to ambiguity; a corresponding 
notation is used for hole states. The superscripts n and p 
of the densities below indicate that neutrons and protons 
orbital transitions have to be calculated independently, 
so they can be later combined to form the transition po- 
tentials. 



1. Elastic from ground state 

The transition density that describes elastic scattering 
from the 0+ RPA ground state is given by: 



^ ' v/2j2 + 1 ^ft.:ia.(P2,Pl) 
Pi>F. P2>F 



hi<F, h2<F 



(15) 



where the prime on the summation sign is a indication 
that only neutron orbitals are to be used for the neutron 
density, and only proton orbitals for the proton density. 
Although in the elastic scattering the occupied orbitals 
do not correspond to holes, we still associate such states 
with the indices hi and /12 in our notation to maintain 
consistency. 

The ground-state density is easily evaluated from the 
above in the Hartree-Fock approximation. The Z coeffi- 
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cients for the particle states vanish, and those for the oc- ^tt(^2, hiY ' — (1)"^^ ^2 + 1 1 



cupied states are simply eZas(^2, hi) = J^i^/n V2j2 + 1- 



Using this in Eq. (|15|) we obtain that the density is 

'(2J/. + 1) 



"(P).9- 



■(r) 



E 



Mr) 

r 



(16) 



^. From ground state to excited states 

After calculating the corresponding spectroscopic am- 
plitudes (Z coefficients) we have that the transition den- 
sity needed for scattering from ground to an excited state 
described by a single RPA boson is 



47r V2 J + 1 

Y: ' [^pV + (-1)^ y,r] ^'^-'^ 



p>F 
h<F 



(17) 



^UmYL{h)^Ss]j\\h\jH)^'^ 

where the initial (ground) state has spin 0. 
The single-particle matrix elements are: 

{i2\mYLy^Ss]j\\k\n)^^{~it 

X hhiihLSJ (^^ o) 1 ^ ^ ^ 



(18) 



The following symmetry relation results from interchang- 
ing initial- and final-state quantum numbers: 

{h\jl\\[YLXSs]j\\h\32)^ 

{-l)'^+'^'^''-'Kh\mYL X Ss]j\\h\ji). (19) 

3. Between excited states 

For scattering involving two excited states, the spectro- 
scopic amplitudes may be divided into two components, 
identified with the superscripts (1) and (2). The parts 
of the spectroscopic amplitudes with superscript (2) are 
present only if the initial and final states are the same and 
J — Q. After a great deal of angular momentum algebra, 
we find the parts of the spectroscopic amplitudes: 



z/,(p2,Pi)(^^ = (-1)'^'-^ v/2j7TT 5](-i; 



32+ jh 



Jf J Ji 
jl jh 32 



(20) 



^P2h 





J 


■'■11 


\j2 


3h 





X 



32 jp jl 



phi 



ph2 



(21) 





J 


■'■11 


[31 


jp 





and 



Z,-^(p2,Pi)(2) = 5 



Nf,JfMf,Ni.JiMi 



Z,{{h2,hi 



l(2) 



X SjM,00 Z^tt%asiP2.Pl) 



— ONfJfMf,NiJiMi 

X SjM,OQ Z^tt%as{h2,hi) 



(22) 



(23) 



where Ni and Nf are the principal quantum numbers of 
the initial and final states, respectively, and Ji and Jf 
are their corresponding angular momenta. 
The corresponding transition density is: 



2 J, + 1 



PlsIW- 4^ y 2j^ + i 



x| Y.' [zit{P2,Pii^'' + Zt,{p2,pif 

I pi>-F, P2>F 

xz'-'^ {l2\32\\[YL{it) X Ss]j\\li\3i) ^ ^ 
+ Y' \zUh2,hi)^'^+Z,i{h2,hif^ 



hi<F, h2<F 



x^''-'' {l2y2\\[YL{rt) X Ss]j\m3i) 
(i)h2{r) 4>hi{r) 1 



(24) 



where the reduced matrix elements are given by Eq. ( 18 1. 
For transitions between different initial and final states, 
only the Z coefficients with superscript (1) appear. When 
the initial and final states are the same, those with super- 
script (2) also appear, but only for zero angular momen- 
tum transfer (J = 0). This case corresponds to elastic 
scattering and we see that the density is the same as for 
elastic scattering from the ground state, but with a cor- 
rection term given by the coefficient with superscript (1). 
When J is nonzero and the initial and final states are the 
same, only the coefficient with superscript (1) appears 
and this can connect different magnetic substates when 
L is even and overall angular momentum is conserved. 
This is a reorientation effect. 



D. Transition Potentials 

The transition potential between any pair of levels of 
a given target nucleus is obtained by a linear combina- 
tion of the transition potentials for protons and neutrons 
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of the target, which are calculated by folding an effec- 
tive nucleon-nucleon interaction with the corresponding 
transition density. This linear combination depends on 
whether the projectile is a neutron or a proton. If v^^{q) 
is the Fourier transform of the effective interaction and 

Pls'/^'^^^^i) Pls/^'''^\i) transition densi- 

ties in configuration space for protons an neutrons re- 
spectively, then the transition potential Uj^gj'^^{q) for a 
(n,n) reaction is 



Tq, fi. 



,Sl\Tq.. fi.ir: 



^Sl^pTq. MP) 



'LSJ 



(25) 



For a (p,p) reaction, the transition potential is 

(26) 

Our scattering effective nucleon-nucleon interaction is 
of Gaussian shape, with parameters matched to the vol- 
ume integral and r.m.s. radius of the M3Y interaction at 
40 MeV; it includes a knock-on exchange correction [35] : 



(27) 



In the momentum space, the explicit form of the effective 
interaction is 



(28) 



where we used the values V§ = -24.1921 MeV and 
/io = 0.7180 fm~^ for the isoscalar part of the interaction 
and = 11.3221 MeV and ^1 = 0.7036 fm^^ for the 
isovector component. We do not include any imaginary 
part in this effective interaction, as our aim is to include 
all nonelastic excitations explicitly in our model. Using 
the Fourier transform v'^{q) of the effective interaction, 
the configuration-space transition potential is 



dqq'jL{qr)v^{q)pll/\q). (29) 



For the reaction calculations coupling to QRPA states, 
we used the single-folded potential constructed by folding 
the interaction of Eq. (28) with the ground-state den- 



sity from the HFB calculation as the bare potential in 
the elastic channel. For simplicity, this was also used 
as diagonal potential for all excited states as well. For 
RPA couplings we used Eqs. (16 1 and (20)-(24| to obtain 



the diagonal potentials. Diagonal potentials describe the 
elastic scattering for the g.s. or an excited state, while 
the off-diagonal potentials provide the couplings between 
the different states of the nuclei. 



nuclear reaction, all possible couplings between channels 
(elastic, inelastic, transfer, etc.) should be explicitly con- 
sidered. However, in practice the number of couplings has 
to be limited, taking into account the relevance of each 
channel. A criterion that may be adopted to define this 
limit is to consider the corresponding order of the cor- 
rection to the elastic scattering. To explore the relative 
importance of the various contributions to the reaction 
cross section, we carried out a series of calculations in- 
cluding three different sets of couplings: A) only inelastic 
couplings from and to the ground state, described in Sec. 
Ill A B) same couplings as in A) but also couplings be- 
tween the excited states, described in Sec. |III B[ and 



C) same as in A) but also couplings to transfer chan- 
nels leading to the formation of a deuteron, described 
in Sec. Ill C[ Table |l] illustrates the different couplings 
that were explicitly included in the three sets of calcula- 
tions (A, B, and C), and the corresponding order of the 
correction from elastic scattering. The contribution of 
the third-order processes included in calculation B) was 
found to be negligible, as shown in Sec. |IIIB| Hence, 
other third-order (i.e. couplings between inelastic and 
transfer channels) and fourth-order (deuteron break-up 
channels, multi-step transfer via continuum) couplings 
were neglected. 



TABLE I: Summary of the channels that were explicitly cou- 
pled in the calculations A, B, and C. The corresponding order 
of the correction relative to elastic scattering is also shown. 



Couplings present 


A 


B 


C 


Correction 


Excited states with g. s. 


Yes 


Yes 


Yes 




order 


Between excited states 


No 


Yes 


No 


3rd 


order 


Transfer with deuteron channel 


No 


No 


Yes 


2nd 


order 



Truncations in channel space are necessary in any prac- 
tical calculation. In our calculations, we depend only on 
real components in the scattering potentials, as we aim 
to reduce the number of phenomenological parameters to 
a minimum. We therefore exclude all channels coupled 
to the elastic at second order and above, and use the 
doorway approximation so that the imaginary parts of 
the optical potentials in the first-order channels hardly 
affect the scattering in the elastic channel. Any imagi- 
nary potentials introduced into our models are only for 
computational convenience. They allow us to temporar- 
ily replace couplings already studied, and to focus on 
additional effects. 



III. COUPLED REACTION CHANNELS 
CALCULATIONS 

According to the coupled reaction channels (CRC) for- 
malism, to correctly account for the effects of the com- 
petition between the different processes occurring in a 



A. Inelastic Coupled Channels 

We initially performed coupled channels calculations 
for reactions involving protons and neutrons scattered 
by the nuclei ^"Ca, ''^Ca, ^^Ni, ^OZr, and ^''''Sm, cou- 
pling the ground state to all levels with excitation energy 
(E*) lying below some limit. Such excited states were ob- 
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tained according to the QRPA model in a box of 15 fm. 
The QRPA states above the particle emission threshold 
are used to approximate exact scattering waves. Re- 
cent studies have shown that such wave functions contain 
large density distributions outside the nuclear radius . 
When used in reaction calculations they accurately rep- 
resent the continuum I37H39| inside the maximum radius. 
Thus, processes containing one nucleon in the continuum 
(plus the inelastically scattered projectile) are included 
in our model. Several applications of this approach can 
be found in the literature for nuclear structure [JD] and 
reactions J?T] problems. Results for the nonelastic ab- 
sorption for neutron-induced reaction on a ^'^Zi target, 
corresponding to the inelastic couplings from g.s., are 
shown in Figure [l] (solid line). 



B. Couplings Between Inelastic Channels 

Couplings between excited states, as predicted by the 
RPA model, were explicitly considered with transition 
densities given by Ec[s. (p0|-p4|) 
by 90 Zr 



200 



for nucleons scattered 
le reaction cross-section 
s^Zr at 



In Figure [Ij we show t 
as a function of partial wave for the reaction n 
scattering energies of 10 MeV (upper panel) and 20 MeV 
(lower panel), where all RPA states lying below 20 MeV 
were coupled. For each i?iab we compare the calcula- 
tions considering only couplings to and from the ground 
state with calculations that also include couplings be- 
tween excited states (with a maximum value, Lmax, for 
the transferred angular momentum between them). It 
is observed that, although for i^iab = 10 MeV there are 
small but noticeable differences between calculations, for 
higher energies the curves are almost undistinguishable. 
Despite the fact that, above 10 MeV, the couplings be- 
tween excited levels do change the cross-sections of the in- 
dividual states (mostly small changes, although for some 
few channels the cross-sections may differ by a factor up 
to 20%), the overall sum among all states remains un- 
changed. This supports the validity of the doorway ap- 
proximation, according to which the total flux leaving the 
elastic channel to all possible first-order channels is in- 
dependent of what happens afterwards: a nucleon later 
might escape as a free nucleon, the flux might equilibrate 
to compound-nuclear resonances, etc. The fact that cou- 
plings between excited states were found to provide neg- 
ligible contribution to the absorption summed over all 
states, for scattering energies above 10 MeV, allowed us 
to disregard them, and other higher-order multi-step pro- 
cesses, in the subsequent CC calculations. 



C. Coupling to Pick-up Channels 

Pick-up channels play an important role in nucleon- 
nucleus scattering |42H44] . We performed coupled reac- 
tion channels (CRC) calculations that included all open 
channels that describe the formation of a deuteron by 



I ' I ' r 

Couplings to/from g. s. only 




Partial Wave 



FIG. 1: (color online) Reaction cross-section as a function of 
the partial wave for the reaction n -I- ^"Zr. Couplings to all 
RPA states below 20 MeV were included. The solid lines rep- 
resent calculations without couplings between excited states. 
The dashed lines correspond to the results when transitions 
between excited states are explicitly calculated, with maxi- 
mum transferred angular momentum Lmax = 2. The dash- 
dotted lines correspond to the same calculation as the dashed 
ones but with L,nax ~ 4. 



picking up the appropriate nucleon from occupied levels 
in the target. For transfers, we approximate the HFB 
target states by bound single-particle states in a Woods- 
Saxon potential, with the radii fitted to reproduce the 
volume radii and Fermi energy obtained by the HFB cal- 
culations. The real diffuseness and spin-orbit parameters 
were taken from Koning-Delaroche optical potentials ^45] 
at i?iab = 0, with spin-orbit radii adjusted by the same 
factor used to fit the volume part to HFB radii. The 
values used are shown on Table HIl 

To overcome numerical limitations, we coupled explic- 
itly only to the transfer channels, incorporating all in- 
elastic effects in an inelastic optical potential obtained 
from coupling only to inelastic channels. The imaginary 
component of this inelastic optical potential corresponds 
to the Koning-Delaroche [45 optical potential renormal- 
ized to account only for the inelastic absorption obtained 



in Calculation A) (described in Section HI A). 

CRC calculations require, in addition to the scatter- 
ing potentials in the incoming channel, a scattering po- 
tential between the deuteron and the remaining target. 
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TABLE II: Parameters of the binding potentials used for the 
coupled reaction channel (CRC) calculations 



Reaction 


a-vifrn) 


Vso (MeV) rsoifm 




'"'Ca(n,d) 


0.67 


5.85 


1.10 


0.59 


4°Ca(p,d) 


0.67 


5.80 


1.12 


0.59 




n fi7 


5.88 


1.13 


u.oy 


«Ca(p,d) 


0.67 


5.82 


1.08 


0.59 


5»Ni(n,d) 


0.67 


5.91 


1.08 


0.59 


5»Ni(p,d) 


0.67 


5.86 


1.10 


0.59 


9''Zr(n,d) 


0.66 


6.01 


1.11 


0.59 


««Zr(p,d) 


0.66 


5.97 


1.07 


0.59 


"''Sm(n,d) 


0.66 


6.19 


1.10 


0.59 


"*Sm(p,d) 


0.66 


6.16 


1.08 


0.59 



We adopted the Johnson-Soper [46] prescription as it in- 
cludes the effects of deuteron breakup in adiabatic (sud- 
den) approximation. In this prescription, the deuteron 
potential is the sum of the individual neutron and pro- 
ton potentials with the target. For the real parts we used 
the diagonal transition potentials of the corresponding 
nucleon-nucleus reaction and, for the imaginary parts, 
the sum of the imaginary parts of the Koning-Delaroche 
pS] optical potential for protons and neutrons. That is, 
fitted parameters are used in the imaginary part of the 
deuteron potential. A phenomenological deuteron poten- 
tial, such as the one proposed by Daehnick et al. [47] . 
was not used since it would introduce fitted parameters 
also in the real part of the deuteron potential. An itera- 
tive method of eliminating this need of phenomenological 
parameters is currently being investigated but we leave 
for future work to calculate the deuteron and nuclear po- 
tentials self-consistently. 



IV. RESULTS 

To assess the success of our large-scale coupled- 
channels approach, we compare the calculated reaction 
cross section to that obtained by the Koning-Delaroche 
optical potential [45], which is one of the best nucleon- 
nucleus phenomenological optical potentials available, 
henceforth referred to as c^^- 

We examined the convergence with respect to maxi- 
mum excitation energy by performing CC calculations 
that explicitly couple to all QRPA excited states below 
the cutoff energies of 10, 20 and 30 MeV, for the incident 
energies = 10, 20, 30 and 40 MeV. We found that 
convergence of the inelastic calculations requires coupling 
of all excited levels below the scattering energy (i.e. all 
open channels) . This behavior is illustrated in Figure [2j 
for neutrons scattered by ^^Ni. As it is seen in Figure 
[2] upper panel, coupling to more highly excited states 
results in larger reaction cross sections. In principle, vir- 
tual couplings between the g. s. and closed channels can 
affect the reaction cross sections by changing the real, 
but not the imaginary, part of the OMP. However, in the 



lower panel of Figure [2] we observe that the couplings to 
states lying above the scattering energy do not noticeably 
contribute to the absorption. 
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FIG. 2: (color online) Reaction cross-section as a function of 
the partial wave for the reaction n -|- ^*Ni at Eiab = 30 MeV 
(a) and 20 MeV (b). 

For reactions having protons instead of neutrons as 
projectile, such inelastic convergence is reached with 
lower cutoffs, as is observed in Figure [3] for p -I- ^^Ni. 
This is due to the fact that a charged particle has to over- 
come the Coulomb barrier, which decreases the amount 
of energy available to excite target states. Thus, states 
close to -Eiab that would correspond to open channels in 
the case of an incoming neutron become effectively closed 
for proton projectiles. This behavior is observed for each 
partial wave as well as at all energies, as is illustrated for 
p -I- ^^Ni in Figures [3] and |4j respectively. Figure [3] shows 
the reaction cross section as a function of the partial wave 
for incident protons of 30 MeV. We observe that the dif- 
ference between the results when coupling to states up 
to 20 and 30 MeV is much smaller than the one for the 
similar system having neutrons as projectile, shown in 
Figure [2] (upper panel). We also observe in Figure |4] that 
the effects of coupling to states lying between 20 and 30 
MeV become noticeable only for E\^\^ > 27 MeV, which 
is consistent with the Coulomb barrier height of « 7.4 
MeV, for this reaction. 

Although the reaction cross section increases with the 
number of coupled inelastic states, to the limit where 
all open channels are coupled, Figure |4] shows that such 
inelastic couplings account only for a small fraction (w 
39% at Ei^b = 30 MeV) of crg'^. However, after including 
couplings to the pickup channels through the CRC calcu- 
lations, a large increase is found, approaching cr§^ and 
the experimental data, as can be seen in Figure [5] An 
even better agreement can be obtained after we include 
the non-orthogonality terms '45', p. 226] in the CRC cal- 
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FIG. 3: (color online) Reaction cross-section as a function of 
the partial wave for the reaction p + ^*Ni at Eiab = 30 MeV. 
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FIG. 4: (color online) Total reaction cross-section as a func- 
tion of the incident energy for the reaction p -I- ^*Ni, for the 
different inelastic calculations. The short-dashed line shows 
the results using the Koning-Delaroche [45j optical potential. 
The lines serve as guides to the eye as calculations were per- 
formed only for ^i^b = 10, 20, 30 and 40 MeV. 



culations, also shown in Figure [Sj Non-orthogonality cor- 
rections arise because at small radii the deuteron bound 
state is not orthogonal to bound states occupied in the 
target. 

In Figure [6] we directly compare our results for p -|- 
^^Ni, ''^Ca with experimental data taken from the liter- 
ature. We show the reaction cross sections obtained by 
coupling only to inelastic open channels (i. e., states with 
E* < -Eiab)j and also by including couplings to pickup 
channels (with non-orthogonality corrections) in addi- 
tion to the inelastic channels. As observed in Figure |6j 
despite the important role of the inelastic channels, the 
most significant contribution to absorption comes from 
the pickup channels. We achieve a good description of 



the experimental data after appropriately including cou- 
plings to all aforementioned processes. 

Within the doorway approximation, we ignore explic- 
itly the absorptive effects of couplings to resonances in 
other channels, such as compound nucleus formation. 
It is thus physically expected that, to compensate for 
the missing absorption from the neglected resonances, 
a damping in the non-elastic channels is necessary. In 
the calculations mentioned above, a damping imaginary 
component, corresponding to the imaginary part of the 
Koning-Delaroche optical potential 45J, was added to 
the inelastic diagonal potentials. Imaginary potential 
components had already been included in the transfer 
channels. For comparison purposes, we show the result 
obtained when this imaginary component in the inelastic 
channels is not added (undamped) for p -I- ^^Ni in the top 
panel of Figure [6j There is noticeable difference between 
damped and undamped calculations when including only 
inelastic couplings (the damped line approximately re- 
produces the lower limit of the undamped oscillations), 
but there is almost no discrepancy when the pickup chan- 
nels are also coupled. The observed differences between 
the damped and undamped calculations indicate that the 
doorway approximation adopted, although very good, is 
not perfect. 

A finite-range interaction in a HFB description of the 
target structure was also considered, as described in Ref. 
|53| . For reactions of nuclcons scattered by ^°Zr, the re- 
action cross section results using the QRPA model with 
the SLy4 functional were found to be practically equiva- 
lent to the results found using RPA states and transitions 
with the Gogny DIS force [SI], as it is illustrated in Fig- 
ure [7] This was observed despite the proton pairing gap 
of 1.2 MeV of ^"^Zr. The fact that, below a given excita- 
tion energy cutoff, there are generally many more QRPA 
states than RPA ones was not significant. In Figure [8] 
we show the total reaction cross sections for p -I- ^°Zr as 
a function of the incident energy obtained by coupling 
to RPA states and using the Gogny DIS force, obtain- 
ing again good agreement with experimental data after 
including the pickup channels. 

This work focuses mostly on reaction cross sections, 
which test the modulus of the S-matrix elements. Ad- 
ditional insights can be gained from elastic angular dis- 
tributions. Preliminary calculations of these give rea- 
sonable agreement with measured cross sections. As an 
example, we show in Figure [9] predictions for the elas- 
tic cross sections of 40 and 65 MeV protons scattered 
by a ^°Zr target. The dash-dotted lines represent cal- 
culations performed within our model, coupling to all 
QRPA inelastic and deuteron channels. A phenomeno- 
logical spin-orbit component from [45] was added, so the 
analysis of our model would be, for now, limited to the 
central components of the OMP. The single-folded po- 
tential obtained from the HFB ground-state density was 
used as the central bare potential, as described previously 
in Section |II D| The solid lines correspond to the same 
calculation but with a modified bare potential instead. 
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FIG. 5: (color online) Reaction cross-section as a function of 
the partial wave for the reactions n, p + *'''''*Ca, ^^Ni, ^*''Sm, 
for the different calculations at Siab = 30 MeV. The results 
shown include couplings to the inelastic states lying below 30 
MeV (dashed green lines), to the inelastic and transfer chan- 
nels (dash-dotted blue lines) and to the inelastic and trans- 
fer channels with non-orthogonality corrections (solid black 
lines). The Koning-Delaroche |45j optical model calculations 
are shown as short-dashed red lines. 



which has a ^-^18% reduced internal depth while main- 
taining the same intensity in the surface region. This 
way, the internal depth of this new modified potential 
approaches the phenomenological one while the reaction 
cross sections remain essentially unchanged. As observed 
in Figure [9j these results are in better agreement with 
experimental data, especially for smaller angles, provid- 
ing an additional improvement from the calculation with 
the unmodified potential. This illustrates how the angu- 
lar distributions are sensitive to the effective interaction. 
This sensitivity can serve as a method to not only evalu- 
ate how realistic the initial structure assumptions were, 
which were consistently extended to the calculation of 
reaction observables, but also to identify important fea- 
tures of the effective interaction adopted. For example, 
the fact that a better agreement with experimental data 
is achieved by reducing the internal depth of the bare po- 
tential indicates that density-dependence effects, which 
have been ignored explicitly, should be relevant for fu- 
ture refinements of the model. Considering the energy 
dependence of the interaction, or implementing an ex- 
plicit treatment of the exchange terms would also im- 
prove the agreement of the model predictions with the 
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FIG. 6: 



(color online) Total reaction cross-section as a func- 

58t 



tion of the incident energy for the reactions p -I- Ni (a) and 
p -I- ''^Ca (b). The results are shown for calculations that 
include couplings to the damped inelastic states lying below 
the scattering energy (dashed green lines), and to the inelas- 
tic and transfer channels with non-orthogonality corrections 
(solid black lines). The Koning-Delaroche optical model 
calculations are shown as short-dashed red lines. For the ^**Ni 
target (a), we also show the results of two calculations that do 
not include damping, one for inelastic couplings (dotted ma- 
genta line) and one for inelastic and pickup channels (dash- 
dotted blue line). Data from Refs. |49H52] . 



experimental data. The result from a phenomenologi- 
cal optical potential (351 calculation is also shown in the 
figure for comparative purposes. 



V. CONCLUSION 

We have calculated the reaction cross-sections for nu- 
cleon induced reactions on nuclei ^°''**Ca, ^^Ni, ^°Zr and 
^^^Sm. This was done by explicitly calculating the cou- 
plings to all relevant transfer and RPA and QRPA inelas- 
tic channels. It was found that this inelastic convergence 
is achieved when all open channels are coupled. Inelastic 
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FIG. 7: (color online) Comparison of the reaction cross- 
section, as a function of the partial wave, between calcula- 
tions using RPA states and transitions with the Gogny DIS 
force [54] (dashed lines) and using the QRPA model with the 
SLy4 functional. The results shown are for neutrons (a) and 
protons (b) scattered by the target ^"Zr, both with Ei^h = 
30 MeV. Couplings to all states below 30 MeV were included, 
according to each model. Transfer channels were not included 
in the comparison. 



couplings account for an important part of the reaction 
cross section, but still most of it corresponds to couplings 
to the deuteron formation pickup channel. The effect of 
couplings between excited levels was studied, leading to 
the conclusion that they did not contribute significantly 
to the nucleon-nucleus reaction cross sections, and thus 
could be ignored. We were able to obtain reaction cross 
sections that were in good agreement with phenomeno- 
logical optical model results and experimental data, by 
means of coupled channels and coupled reaction channels 
calculations, with non-orthogonality corrections. Ab- 
sorptive effects of resonances from other channels appear 
to be negligible when the one-nucleon pickup channels are 
coupled in addition to the inelastic channels. Preliminary 
calculations of elastic angular distributions were also in 
good agreement with experimental data, for small an- 
gles. This demonstrates the applicability of the doorway 
approximation, made evident by the negligible difference 
in absorption observed between calculations considering 
damped and undamped inelastic couplings, after also in- 
cluding the pickup channels. Future work on couplings 
between different types of nonelastic processes will in- 
clude higher-order corrections. 
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Appendix A: Fourier-Bessel Expansions 

The following discussion is based on the treatment of 
Petrovich [57j, which has been applied to inelastic scat- 
tering by Petrovich, Carr, and McManus in Ref. [55] . 

We consider the Fourier-Bessel series expansion of a 
function that can be expressed in the form 



/(r) = fjM{r)Y:j,,{v), 



(Al) 



which is defined within a radius R. We refer to the radial 
part /,/m(^) as a partial- wave function. In applications 
the function /jm is frequently independent of M and is 
then labeled /j; we retain the more general form. The 
Fourier transform is defined as 

/(q) ^ j dv e±^i V(r), (A2) 

where the domain of the radial part of the integral is to 
R. Petrovich employs the lower (— ) sign; we retain both 
possibilities. 

By using the expansion of the exponential term of Eq. 



(A2) 



e±^q-"- = Ai,Y.{±zy jj{qr) r,V(f )r,M(q), (A3) 



JM 



and the definition of the dot produc t of spherical tensors 
(SQ) and inserting the result in Eq. (A2) we find: 



/(q) = (±i)^/jM(g)r;M(q), 



(A4) 
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FIG. 9: (color online) Angular distribution for the elastic 
cross section for the reaction p + ^'^Zr at scattering energies 
of 40 (a) and 65 MeV (b). The dash-dotted and solid lines 
correspond to the predictions of our model using the poten- 
tial from HFB densities and from a slightly modified bare 
potential (see text), respectively. The Koning-Delaroche [45) 
optical model calculation is shown as the short-dashed line. 
Data from Refs. i55. .56. . 



/jM(g) =47r/ dr r'jj{qr)fjMir), (A5) 
Jo 



which is the partial-wave form of the transformation; we 
refer to the function /jAf('z) as the partial- wave function 
in momentum space. 

In the hniit i? — !■ cx) it is easy to see that the expres- 
sion for the reverse transformation is the usual Fourier 
integral; that is, 



fir) 



(27r)3 



rfq e^"i V(q), 



(A6) 



where the integral extends over all momentum space, or 
in partial-wave form 



fjAiir) 



477 



dq q^ij{qr)fjM{q). 



(A7) 



(27r)-^o 

If the partial wave expansion of the function to be 
transformed is defined in terms of Yjm instead of Y 
the above discussion is unaltered 



That is, if Ff^ 



is re- 



placed by Yjm in Eq. (Al), all expressions are identical 



except for the same replacement in Eq. (A4) 
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